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ABSTRACT 

We use published measurements of the transmission power spectrum of the Lyman a forest to constrain 
several parameters that describe cosmology and thermal properties of the intergalactic medium (IGM) . 
A 6 parameter grid is constructed using Particle-Mesh dark matter simulations together with scaling 
relations to make predictions for the gas properties. We fit for all parameters simultaneously and identify 
several degeneracies. We find that the temperature of the IGM can be well determined from the fall-off 
of the power spectrum at small scales. We find a temperature around 2 x 10'' K, dependent on the 
slope of the gas equation of state. We see no evidence for evolution in the IGM temperature. We place 
constraints on the amplitude of the dark matter fluctuations. However, contrary to previous results, the 
slope of the dark matter power spectrum is poorly constrained. This is due to uncertainty in the effective 
Jeans smoothing scale, which depends on the temperature as well as the thermal history of the gas. 

Subject headings: cosmic microwave background — methods: data analysis 



1. INTRODUCTION 



shifts around z ~ 3, which can place strong constraints 



Th a past decade has seen remarkable progress in our 



under rtnnding of the L; 



foroat. Compmioon of the 



on cosmology and the na.ture of dark matter ( Croft et alj 
19981, ICroft et al. 1999| , [White fc Croft 2000| , P?arayanan| 



absor ption lino data and numerical simulations hae load to 
a clear picture for the forest sometimes called the fluctu- 

[Hernquist 



et al. 2000[ ). Croft et al. (1998) pioneered the approach of 



ating Gunn Peterson effect (Cen et al. 1994 



> . l et al. 1995|, [Zhang et al. 1995|, [Mhalda-Escude et al 



19961, IMuecket et al. 1996| , [Wadsley fc Bond 1996| , [Ihe- 



uns et al. 1998). In this picture most of the absorption 
is produced by low density unshocked gas in the voids or 
mildly overdense regions in the universe. This gas is in 
ionization equilibrium and traces broadly the distribution 
of the dark matter, but is also sensitive to its equation of 
state. Simple semi-analytic models based on these ideas 



ing the main features of numerical simulations (|Bi et al 


1992 




Reiscneggcr & Miralda-Escude 1995, 


Bi & Davidser 


1997 




Gnedin fc Hui 1996, Croft et al. 1997 


p3ui fc GnediE 


1997a 


, Hui et al. 1997bD. 



The main ingredient that determines the absorption in 
the forest is the distribution of the dark matter so the 
forest can be a very powerful probe of cosmology. The 
probability distribution of the transmission has also been 



computed and successfully compared to the data ( Ranch et 
al. 1997| , ^sser fc Haehneh 2000i ^cDonald et al. 19991) 



It has been shown that the moments of this distribution 
can be predicted analytically using the kno wn scalings for 



the matter and si mple ideas of local biasing (Gaztanaga & 



Croft et al. 1999) 



inverting the shape of the mass power spectrum directly 
from the 1-D power spectrum of the flux or transmission. 
The amplitude of the mass power spectrum is then ob- 
tained by comparing with simulations. Hui (1999) pointed 
out the possibility of a bias in the recovered shape due 
to redshift distortio ns, and suggested a modified inversio n 
technique (see also McDonald & Miralda-Escude 1999). 
Most work so far focused on cosmological information that 
could be obtained from the large scale forest power spec- 
trum. 

In this paper we obtain cosmological information from 
the transmission power spectrum of the forest using a dif- 
ferent approach. Rather than trying to extract the power 
spectrum of the dark matter by applying some inversion 
technique to the data, we take the observed transmission 
power spectrum as is, and simply compare it with predic- 
tions from a wide range of models. This is in part the 
approach taken by McDonald et al. (1999), but the type 
of models they examined had a fixed thermal evolution. 
Here, we construct a grid of models described by 6 param- 
eters and constrain all parameters simultaneously using 
a likelihood analysis analogous to what has been devel- 
oped for cosmic microwave b ackground (CMB) data (eg. 
Tegmark fc Zaldarriaga 200C| ). 



Perhaps the most important application of the forest is 
to measure the power spectrum of the dark matter at red- 



Our analysis allows us to use the information on smaller 
length scales. In particular, the sharp drop in thepower 
spectrum of the forest on small scales (see figure |l]) is a 
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Fig. 1. — Transmission power spectrum from from McDonald et 
al. (1999) together with our best fits. 



direct consequence of thermal broadening. The fact that 
the small scale power could be se nsitive to the temper ature 
has been noted in the past (eg. rhcuns ct al. 2000 ). We 
are able to place constraints on the temperature of the 
IGM at high redshift. 

There has been some controversy in the literature as 
to what the equation of state of the intergalactic medium 
(IGM) is as determined by the width of absorption lines. 
Different studies have reached different conclusions about 
the time evolution of the gas temperature and its eq uation 
of state dSchaye et al. 1999| , [McDonald et al. 200o|) . Our 
results can be directly compared to those obtained by these 
other techniques. 

The outline of the remainder of the paper is as follows: 
we describe the data used in our analysis in section ^, de- 
scribe our method in ^ present constraints on parameters 
in ^, suggest ways to improve our analysis in the future in 
H and conclude in ^. 



2. DATA 

The most detailed measurements of the transmission 
power spectrum of the forest in the literature have been 
presented by McDonald et al (1999). We will use these 
data for our analysis. Figure |l| shows the data together 
with the best fit models in our grid. The data points and 
error bars were taken straight from the tables in McDon- 
ald et al. (1999), except for the last two points in the plot 
which were only used as upper limits because of concerns 
about metal line contamination. 

In addition, we used the mean transmission as deter- 
mined by McDonald et al. (1999) with their quoted errors 
(table I). 

3. METHOD 
3.1. Model Predictions 



z 


(F) 


4 


3.89 
3.00 
2.41 


0.475 ±0.021 
0.684 ±0.023 
0.818 ±0.012 


0.1293 ±0.0030 
0.1174 ± 0.0056 
0.0789 ±0.0068 



Mean {F - 
F = e"" 



Table 1 

= (f)) and varia nce ( ap) of the transm ission 
as measured by mcdonald et al. 1999 for 



SEVERAL REDSHIFT BINS. 



The model predictions were done using the PM m odel 
for the forest ( proft et al. 1998^ [Meiksin fc White 2000| ) . In 
this model PM simulations are run to compute dark matter 
densities and velocities. The gas density and temperature 
are calculated using simple scaling relations inspired by the 
results of full hydrodynamics simulations. Our objective 
is to constrain the parameters of these scaling relations. 

We ran 7 PM simulations using 128^ particles with a box 
size 16/i~^Mpc (corresponding to 3200km/sec at redshift 
z = 3) of the standard Cold Dark Matter (SCDM) model 
with different spectral indices n = (0.4,0.7,0.9, 1, 1.1, 1.3, 
1.7) for the power spectrum of initial density perturba- 
tions. Having a non-vanishing cosmological constant will 
not significantly alter our results because fi™ is close to 
1 at the relevant redshifts regardless of whether A = - 
the main determining factor is instead the power spec- 
trum at certain scales, which we will constrain. We 
stored the density and velocity at expansion factors a = 
(0.05, 0.065, 0.085, 0.11, 0.14, 0.19, 0.24, 0.31, 0.41). The ex- 
pansion factors are used to span a range of power spectrum 
normalizations . 

The baryon density is obtained by smoothing the dark 
matter density; the smoothing mimics the effect of pres- 
sure forces. Following Gnedin & Hui (1998) we adopt 
(5b(k) = Wf{k)ScDM{^) where W{k) is a Fourier space 
smoothing, which we take to be Wf{k) = exp(-(fc/fc/)2). 
As discussed in Gnedin & Hui (1998), kf is proportional 
to the Jeans scale but the constant of proportionality de- 
pends on the details of the reionization history of the uni- 
verse (see also Nusser 2000). For this reason, we treat kf 
as a free parameter. At redshift z ~ 3, the filtering is 
expected to be around kf ~ 35 hMpc~^ (Gnedin & Hui 
1998). 

The optical depth for each gas element is related to the 
overdensity r using a power law (Hui & Gnedin 1997a), 



a-o 



(1) 



where A = p/p denotes the gas overdensity. The transmis- 
sion is e""^. The constant oq is fixed so that the generated 
spectra have a particular value of the mean transmission 
(which we will call /). We make a grid of mean transmis- 
sions but use the observed measurements (in table |^) as 
a prior in our likelihood. Our grid of mean transmissions 
serves effectively as a grid of ao's. 

For each grid point in the box a temperature is assigned 
to the gas using a simple power law equation of state for 
the gas, 

T = roA". (2) 

We treat Tq and a as free parameters. The power law index 
P in equation (□) is related to a, /3 = 2 — 0.7a. Smoothing 
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due to thermal motions is included when the spectra are 
generated. The absorption produced by each fluid element 
is distributed in velocity space as exp(— (As/6)^)/&-y7r, 
where As is the velocity space separation and the b param- 
eter is given hy b = yj2kT jm^ k, 13 km s^^ for T = 10"* 
K. 

Our parameter vector p — {a,n, kf ,To, a, f) has 6 di- 
mensions. We created a grid of model predictions for each 
choice of parameters in a grid, 

• a = (.05, .065, .085, .11, .14, .19, .24, .31, .41) 

• n = (0.4,0.7,0.9,1,1.1,1.3,1.7) 



• k 



f 



(5, 10, 20, 30, 40, 45, 50, 55, 60, 70, 



• To = (150,200,250,300,350,400,450,500,550,600, 
700) 

• a = (0.0.0.1, 0.2, 0.3, 0.4, 0.5, 0.6) 

• f ^ (.4, .45, .475, .5, .525, .55, .575, .6, .625, .65, .665, 
.684, .7, .725, .75, .78, .8, .82, .85, .9) 

where by convention the scale factor is a = 1/(1 + z), kf is 
measured in hMpc~^ and Tq in (kms~^)^. Thus we have 
a total of9x 7x llx llx 7x 20 = 1067220 combi- 
nations of model parameters. For each point in our grid 
we extract 500 lines of sight randomly from the simula- 
tions, generate spectra and measure the power spectra of 
the forest. We measure the power spectra of the relative 
transmission fluctuations, SF/F where F = exp(— r). 

3.2. Likelihood calculation 

For each model in our grid, we compute the likelihood 
by comparing the power spectra computed form the sim- 
ulations to the points in McDonald et al. (1999). In this 
paper, we will stick to a crude Gaussian approximation, 



£(d; p) cx exp 



1 f d,-Pf,{p) 



(3) 



where i runs over the different data points di, the model 
prediction for that wavevector are P/i(p) and Ui are the 
error bars on each point. By comparing different realiza- 
tions we computed the statistical errors in our model pre- 
dictions, which were approximately 5%. We added those 
in quadrature to the observational errors to compute (7^. 
The addition of this extra source of variance makes very 
little difference because the observational errors are sig- 
nificantly larger. Most certainly systematic errors asso- 
ciated with our simplified model of the forest will domi- 
nate the error budget. We also used the last two points 
(fc > O.lkm^^s) as upper limits. To do so in practice we 
increased cr for those points to cr = 0.3P/, roughly a factor 
of three increase on the quoted error bars and used it as a 
one sided error bar. For models that had more power than 
Pf we computed with the increased error, for models 
with less power, we set the contribution to from that 
point to zero. 

The data points reported in McDonald et al. (1999) 
are for the power spectrum of the transmission SF, and 
this is what we directly compare our predictions with (see, 
however, §H for virtues of using power spectrum of 6F/F 



instead). We furthermore add to_the likelihood a term to 
account for the measurement of F. We rnultiply the likeli- 
hood in ^ by expj— (/ — F)'^/ap) where / is the parameter 
in our grid and F is the observed value at the appropriate 
rcdshift. The error in the mean transmission is denoted 

The full likelihood function is £ = e^^ where is 
simply the chi-squared goodness of fit of the model to the 
data. We have chosen to keep things this simple because 
we are in any case unable to eliminate a major source 
of inaccuracy: most probably there are correlations be- 
tween the estimates of the power on different scales which 
are hard to estimate. We are using 18 bins for the power 
spectrum, thus estimating a covariance matrix requires es- 
timating 153 numbers. Naive estimates of the covariance 
matrix based on a limited number of lines of sight are 
very noisy which translates into artificially low probabil- 
ities for some models, thus they cannot be used in this 
analysis. More theoretical work needs to be done in order 
to address this problem in a satisfactory manner. Perhaps 
one could look for guidance in simulations to build a sim- 
ple parametrised model of the covariance matrix and then 
use the scatter between the lines of sight to fix those pa- 
rameters. Moreover the window function needed to relate 
observed and theoretical power spectra are not available 
for this observation. We choose to use flat windows inside 
each bin but also tried Gaussian windows, which made 
only minor differences. There is ample room for improve- 
ment in this part of the calculation. 

Once we have computed the likelihood for every model 
in our grid we marginalize along one direction at a time 
until we get one- or two-dimensional constraints on pa- 
rameters. We use the method developed in Tegmark & 
Zaldarriaga (2000) for this purpose. 

4. RESULTS 

In this section we report the constraints we have ob- 
tained. We first discuss our results on the equation of 
state of the IGM, then on the relation between baryons 
and dark matter and finally we summarize our constraints 
on the power spectrum of primordial fiuctuations. 

4.1. The temperature of the IGM 

The shape of the power spectrum of the forest has a 
sharp cut-off around k = 0.02 km^^ s which is due to 
the thermal broadening of the lines. Figure ^ shows a se- 
quence of models with increasing IGM temperature. As 
the temperature increases the smoothing becomes more 
important and the power on small scales is reduced. At 
k 0.04 km"-'^ sec the power is down a factor of 2 from the 
extrapolation of the larger wavelengths. This corresponds 
to a temperature T ~ 1/fc '--^ 2 x W^K. Our likelihood 
analysis takes advantage of this dependence to find con- 
straints on the temperature. 

Figure summarizes the constraints we obtained. It 
shows 95% confidence regions]^ in the Tq — a plane. Only 
a's between and 0.6 are examined, because realistic a 



are expected to fall in this range (Hui & Gnedin 1997a). 
We find a degeneracy in this plane, which is caused by the 

^The lines that enclose our allowed regions are actually lines of 
constant x^- particular we use Ax^ = 6.18, which would corre- 
spond to 95% of the area under a 2 dimensional Gaussian. 
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k [sec km"'] 

Fig. 2. — Transmission power spectrum at z = 3 for mod- 
els with varying Tq (Tq = (150, 250, • • ■ , 650)(km/s)^ increasing 
from top to bottom). All other parameters are left constant, 
p = (.24, .7, 50, To, .2, .7). The data points are from McDonald et 
al. (1999). 
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Fig. 3. — Constraints in the T — a plane (95 % as described in 
the text). The top panel corresponds to the mean density while the 
bottom panel is for an overdensity A*. The ellipses with lines at 
45° angle on the top are the equivalent constraints from Schaye et 
al. (2000). and on the bottom from McDonald et al. (1999) 



fact that our constraint is most tight for the temperature 
at a density different from the mean density. Under the as- 
sumption that the temperature density relation is a power 
law, the temperature at any given overdensity A* can be 
calculated from (To, a) as: 

-ToA?. (4) 

Our method to determine the temperature is most sensi- 
tive to T*, so models with the same are all good fits 
to the data. This set of models correspond to lines in the 
To — a plane. 

To investigate what A, is, we run simulations in which 
the temperature-density relation has a step at a fixed over- 
density A jump, 

T = To (if A < Ajump) 

= 2To (if A > Ajump) (5) 

with To = 1.2 lO'^K. In figure ^ we show the change in 
the transmission power spectrum as we change the value 
of Aju7np for a model that fits the observations at redshift 
z = 3. For low values there is no change. When Aj^mp ~ 1 
changes become noticeable but once A jump > 1-8 fur- 
ther increases in Ajump no longer change the power spec- 
trum. At redshift z — 3 the transmission power spectrum 
is therefore sensitive to 0.9 < A < 1.8. Given that our 
method is most sensitive to A's larger that one, the con- 
tours of equal likelihood are tilted in the (To, a) plane in 
figure ^. On the other hand, because the method is in- 
trinsically sensitive to a range of overdensities, it should 
be possible to extract information about a once the errors 
in the measured power decrease. The bottom panel of fig- 
ure H shows our constraints at A, = 1.4 constructed by 
reparametrizing our model grid using equation ^. 
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Fig. 4. — Differences in the power spectrum of the flux of the forest 
at redshift 2 = 3 for models with a step like temperature density 
relation as described in the test. Labels on the right indicate the 
value of A jump 

A complementary method to determine the temperature 
was proposed in Schaye et al. (1999) and implemented in 
McDonald et al. (2000), Ricotti et al. (2000)and Schaye et 
al. (2000). The method is based on the observation that 
the width of the lines in the spectrum has several contribu- 
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tions, from thermal broadening, peculiar velocities and th e 
intrinsic sizes of the structures ( Hui and Rutledgc 1999D . 
The minimum width of the lines is thus a measure of the 
thermal broadening. In the three papers cited above, the 
distribution of widths of lines was used to infer constraints 
on the temperature, two of which are shown as ellipses in 
our plot. We emphasize that those constraints were ob- 
tained fitting lines and determining a minimum width in 
their distribution. Our constraint did not require the fit- 
ting of lines. Although the analyses of Schaye et al. (1999) 
and McDonald et al. (2000) are based on the same under- 
lying idea and very similar data sets, they differ in several 
technical details and agree neither in their conclusions nor 
in the size of their error bars. 

Our results are in good agreement with those of McDon- 
ald et al. (2000) and we have comparable error bars. Our 
temperature is higher than that of Schaye et al. (2000). 
This could signal a bias in one of the methods, an under- 
estimate of the errors or some unaccounted for physical 
effect. For example, to interpret the line fitting results 
one needs to know what gas density is responsible for the 
majority of the absorption. Simulations are usually used 
for this purpose and errors in this determination can lead 
to changes in the inferred temperatures. If the difference 
is not systematic, then a particularly interesting possibil- 
ity is an inhomogeneous equation of state, which could be 
produced by a recent reionization of Helium II (see e.g. 
Madau 2000 for a review). In that case, the temperature 
measured from the cut-off in the line width distribution 
would tend to be lower because the method tends to use 
lines in places with a smaller mean temperature. 

4.2. The relation between dark matter and baryons 

The baryons experience pressure forces, so their distri- 
bution on small scales may not be the same as that of 
the dark matter. In our technique, we model this effect by 
smoothing the dark matter distribution to get the gas den- 
sity using a 3D filter with a constant filtering scale across 
the box (kf). In analytic models, kf depends on the tem- 
perature of the IGM and on its reionization history. The 
smoothing scale is one of the parameters in the grid so we 
marginalize over it when obtaining constraints on other 
parameters. 

Even though we have A:^- as a free parameter to marginal- 
ize over, we still get constraints on the temperature of the 
IGM. It is interesting to understand why that is the case, 
because both effects are a form of smoothing. The key is 
that the kf smoothing of the density field is done in 3D. 
After smoothing, a nonlinear transformation is applied to 
the density to obtain the flux. This transformation shifts 
power between scales, distorting the shape of the trans- 
mission power spectrum. Figure ^ illustrates the effect. 
For large smoothing scales (small kf), the shape of the 
transmission power spectrum is totally wrong. Only when 
the smoothing by kf is subdominant to smoothing pro- 
duced by thermal broadening are the fits acceptable. Note 
that some of our higher kf^s probably approach the res- 
olution limit of our simulations. We obtain lower limits 
on kf, kf > 10,30,25 hMpc"! for z = 2.4,3,3.9 which 
can be compared with the Nyquist cut-off of the simula- 
tion kj^yq — 50hMpc~^. Thus some care is needed when 
interpreting our lower limits, as they may be somewhat 
sensitive to our resolution. In any case, the conclusion 
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Fig. 5. — Transmission power spectrum at z = 3 for models with 
varying kf {kf = (5, 10, 20, 30, 40)hMpc-l from bottom to top). All 
other parameters are left constant, p = (.24, .7, kf, 400, .2, .7) . The 
data points are from McDonald et al. (1999). 



that the 3D smoothing should be subdominant compared 
with what is produced by the temperature should be ro- 
bust against increase in numerical resolution. 

4.3. The power spectrum of mass fluctuations 

Finally we want to consider constraints on the power 
spectrum of mass fluctuations. Figure ^ shows our con- 
straints in the amplitude/spectral index plane. The top 
panel shows directly the contours in our primary variables 
{a,n), the expansion factor of our simulation and the pri- 
mordial spectral index. 

In the bottom panel of figure ^, we have changed our 
contours to the (A^,ne//) plane, where A is the ampli- 
tude and Ueff is the spectral index of the linear power 
spectrum at fc = 0.008 km~^ sec. We did this to compare 
with the constraints from Croft et al. (1999) shown with a 
red ellipse. Although we agree in the range of allowed am- 
plitudes, we disagree in the ability of the data to constrain 
the spectral index. The reason we cannot constrain the 
spectral index from the shape of the transmission power 
spectrum can be traced to our having fc^ as a free param- 
eter. Figure || shows that changes in kf modify the low fc 
slope of the transmission power spectrum. 

In figure ^ we plot two models with very similar trans- 
mission power spectrum but with n = 0.4 and n — 1.6. 
The difference in slope, which should otherwise be notice- 
able at low fc, is compensated by a change in fc/. This 
change in fcy should produce changes in the power spec- 
trum at high fc, which is offset by a change in Tq. Thus the 
power spectrum data alone cannot accurately constrain 
the shape of the matter power spectrum unless other in- 
formation is used to constrain kf. 
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5. IMPROVEMENTS FOR THE FUTURE 
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Fig. 6. — The Constraints on amplitude and spectral index of 
the power spectrum of initial fluctuations. In the top we show con- 
straints in our grid variables. On the bottom we have transformed 
this constraints to amplitude and spectral index of the initial power 
gpcrfviim at h — n nn« km~^ sec to compare with the constraints 
form Croft ct al. 199S| shown with an ellipse (lines at 45° angle). 
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Fig. 7. — Two models with very different spectral index that have 
similar flux power spectra together with the z = 2.4 data. One 
model has p = (.41, .4, 70, 350, .3, .82) with = H and the other 
has p = (.14, 1.6, 20, 250, .2, .82) with = 10. 



The objective of this work was to show that detailed 
analysis of the transmission power spectrum can be used 
to extract useful information about the cosmological model 
and the properties of the IGM at high redshift. In order 
to make progress and extract more information from the 
data, there are several aspects of our work that could be 
improved in the future. 

From an observational perspective, we need to improve 
the determination of the transmission power spectrum. 
Important ingredients will be to obtain a reliable covari- 
ance matrix for the different band powers and to determine 
the window functions that relate the observed power with 
the theoretical predictions. It would also be desirable to 
measure the power spectrum of 5F / F after trend-removal 
rather than the power spectrum of 5F after continuum- 
fitting (Hui et al. 2000), to minimize systematic errors 
due to continuum placement. 

From a theoretical perspective, one should study in more 
detail the relation between baryons and dark matter. Hy- 
drodyna mic simulations as we ll as the hydro-PM approx- 
imation ( [Gnedin fc Hui 1998 ) should be used to under- 
stand and quantify the possible biases introduced when 
extracting different parameters from PM simulations. Our 
method requires exploring a large number of models, so it 
will rely on dark matter only simulations (or the hydro-PM 
variant) for the foreseeable future. 

A very interesting possibility is to try to determine if 
there are spatial fluctuations in the equation of state of the 
IGM. This will require measuring the transmission power 
spectrum and inferring a temperature for individual por- 
tions of spectra along different lines of sight and comparing 
results. A comparison will only be meaningful if one un- 
derstands the expected distribution for these band-powers 
from simulations. It is important to know what one ex- 
pects from cosmic variance alone and how different band 
powers are correlated in models with uniform equation of 
state. 

Another direction which should be pursued is the com- 
bination of the power spectrum analysis with other infor- 
mation, in particular the probability distribution function 
(PDF) of the transmission. It will be interesting to know 
if the PDF can break some of the degeneracies we en- 
countered. It is expected however that some degeneracies 
between parameters will remain (Theuns et al. 2000). 

6. CONCLUSIONS 

We have implemented a simple likelihood method to 
constrain the parameters of a Lyman a forest model from 
transmission power spectrum data. We used a 6 parame- 
ter grid and a marginalization method based on previous 
work for CMB data Tegmark & Zaldarriaga (2000). It is 
important that we fit for the 6 parameters simultaneously, 
as we find several parameter degeneracies that prevent us 
from extracting tighter constraints. 

Perhaps our most interesting constraint is that on the 
equation of state of the gas. Our results are in good agree- 
ment with those of McDonald et al. (2000). We find 
a temperature that is somewhat higher than might have 
been expected and see little evidence for strong evolution 
of the temperature with redshift. Our determinations of 
the mean temperature and of temperature-density slope 
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are highly degenerate. 

As for the constraints on the mass power spectrum, 
we find amplitude constraints that are in agreement with 
those found by Croft et al. (1999). On the other hand, we 
are unable to place accurate constraints on the spectral 
index of the power spectrum, mainly due to an uncer- 
tainty in the relation between baryons and dark matter, 
parametrized by a smoothing scale kf. This uncertainty is 
in part a consequence of our using dark matter only sim- 
ulations to model the forest. Thus if one could construct 
a grid using hydrodynamic simulations, perhaps using the 
hydro-PM method of Gnedin & Hui (1998), some of this 
uncertainty could be removed. However, that part of the 
uncertainty in fcj which is due to our ignorance of the 
reionization history will remain. 

An important finding regarding smoothing is that the 
3D smoothing parametrized hy kf is less important than 
the ID smoothing due to thermal broadening in determin- 
ing the small scale transmission power spectrum. This is 
what allows us to determine the temperature of the IGM 
from the rapid fall-off in the small scale power. A kf small 
enough to affect the fall-off would also modify the shape 
of the transmission power spectrum sufficiently to make it 
an unacceptable match to observations. 
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